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ABSTRACT 

A model for gas outflows is proposed which simultaneously explains the cor- 
relations between the (i) equivalent width of low ionization and Lya lines, (ii) 
outflow velocity, and (iii) star formation rate observed in Lyman Break Galax- 
ies (LBGs). Our interpretation implies that LBGs host short-lived (30 ± 5 
Myr) starburst episodes observed at different evolutionary phases. Initially, 
the starburst powers a hot wind bound by a denser cold shell, which after sa 5 
Myr becomes dynamically unstable and fragments; afterwards the fragment 
evolution is approximately ballistic while the hot bubble continues to expand. 
As the fragments are gravitationally decelerated, their screening ability of the 
starlight decreases as the UV starburst luminosity progressively dims. LBG 
observations sample all these evolutionary phases. Finally, the fragments fall 
back onto the galaxy after ss 60 Myr. This phase cannot be easily probed as it 
occurs when the starburst UV luminosity has already largely faded; however, 
galaxies dimmer in the UV than LBGs should show infalling gas. 
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1 INTRODUCTION 

Recent observations (Adelberger et al. 2003, Shapley et al. 2003 (S03), Steidel et al. 2003, 
Adelberger et al. 2005) based on a combination of quasar absorption-line and faint-galaxy 
techniques applied to the same cosmic volumes have revealed that UV-selected, star form- 
ing galaxies at redshift z ~ 3 (in brief, Lyman Break Galaxies, LBGs) show spectral ab- 
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sorption lines due to heavy elements. These absorption lines, which can be very strong 
(New ^ 10 14 cm~ 2 ), are often called "interstellar" because they resemble the absorption 
lines produced by interstellar material in local galaxies. Nevertheless, these lines appear 
(without exception) to be blueshifted with respect to the galaxy systemic velocity deduced 
from nebular lines. The typical velocity difference is Av ps 250 km s" 1 , but a non-negligible 
fraction of LBGs show differences in excess of 300 km s _1 (Pettini et al. 2002). These absorb- 
ing systems are often interpreted as arising in cool and/or neutral structures (maybe a shell) 
part of an outflow powered by supernovae associated with the detected galaxy star forma- 
tion activity. An additional element in support of this hypothesis comes from the detection 
of redshifted Lya photons backscattered from gas behind the stars. Quite interestingly, this 
emission could trace the receding part of the outflowing shell. Finally, strong Lya absorp- 
tion is produced by the intergalactic gas within \hr x comoving Mpc of most LBGs, but for 
about 1/3 the absorption is weak or absent. The transparency of the medium around some 
LBGs could arise from the gas collisional ionization/heating by the wind (Croft et al. 2002, 
Adelberger et al. 2003, Kollmeier et al. 2003, Bruscoli et al. 2003, Bertone 2005, Kollmeier 
et al. 2006), from a galaxy proximity (photoionization) effect, as pointed out by Maselli et 
al. (2004), or simply reflect cosmic variance and artifacts of the data reduction (Desjacques 
et al. 2004, Desjacques, Haehnelt & Nusser 2005). 

In general, one would expect to observe also infalling gas in the surroundings of LBGs 
arising either from gas accretion along the cosmic web filaments or because at least part of 
the outflowing gas could rain back onto the galaxy if its velocity does not exceed the escape 
speed of the system (a process known as galactic fountain). In fact, the infall ram pressure 
might be of primary importance in confining and reverting the wind expansion (Fujita et 
al. 2004). This situation is clearly appreciated from the inspection of the results of recent, 
high resolution simulation of LBG outflows (Ferrara, Scannapieco & Bergeron; see Fig. 2 of 
that paper). However, observations so far have not found any experimental evidence of such 
physically-based expectation. 

The central aim of this paper is to understand why outflows are routinely observed in 
LBGs whereas gas infall seem to escape detection. This puzzle is possibly related to the 
following scenario. During the first phase of a starburst the galactic wind produces a cavity 
filled with hot gas, enclosed by a cold and dense shell. Afterwards the shell fragments and 
the velocity of the cold debris starts decreasing while the hot gas escapes from the galaxy. 
Whether or not they survive for some time, the UV light from the galaxy should become 
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progressively less absorbed by the debris since their covering factor decreases with time as 
they travel to larger distances from the galaxy. Observations may show a time evolution of 
the equivalent width (EW) of metal absorption lines and dust reddening produced by the 
decrease of the screening effect of the cold debris on the UV light from the galaxy. If the 
debris are not disrupted or ejected from the galaxy we should also observe the gas return 
(infall) process. The relevant question to pose then seems: Can we observe galaxies during 
all these different evolutionary stages? 

A large statistical galactic sample (more than 800 galaxies at intermediate redshifts) 
that can be used to test this model has been compiled by S03. From the analysis of this 
data set a very interesting number of relationships emerge (some of which already noticed 
by S03 and others obtained by correlating their data in a different way), which are listed in 
the following: 

(1) The wind velocity, v w , is proportional to the star formation rate (SFR): v w oc M± 

(2) The Lya emission line EW is inversely proportional to the wind velocity: EWi sa <xl/v w 

(3) The EW of metal absorption lines is proportional to the wind velocity: EW met oc v w 

(4) The reddening is proportional to the wind velocity: E(B — V) oc v w 

Should we interpret these properties as produced by intrinsic variations of the physical 
properties (eg., SFR, mass, size) of LBGs? Or do they rather reflect the different ages of a 
single (broadly speaking) galactic population with uniform properties? The relationships (3) 
and (4) naturally descend from the properties of a time evolving uniform population: as the 
cold debris move outward they are slowed down by gravity /energy losses, their covering factor 
decreases, and the EW and reddening also decrease. Property (2) can easily be understood 
in the same framework. Contrary to the metal absorption lines, the EW of the Lya emission 
line increases with decreasing wind velocity (ie., increasing time from the starburst) since 
the scattering by the outflowing neutral gas is reduced and a larger fraction of the blue wing 
of the Lya line is transmitted. 

Naively, relationship (1) can be interpreted as due to a variation of the energy input 
of the wind for galaxies with different star formation rate. However, a closer look at that 
relationship (Figure [TJ shows that the dependence of v w on the SFR is not as the one 
expected under this assumption. Assuming, for simplicity, a constant density profile for the 
gas in the galaxy core, the velocity is a very weak function of the SFR: v oc Ml with q — 1/5 
(shown by the solid line). Even assuming a more general gas density profile of the galaxy 
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Figure 1. Wind velocity as a function of the LBG star formation rate. The points with errorbars show the data from S03; the 
dashed line is the best fit to the data with a power law v w <x M^ 3 . A pure luminosity effect would produce a much flatter 
power law: v w oc M* / ' 5 (solid line). 

core, p g oc r~^, would produce q = 1/(5-/5) that is too flat to fit the observations for 
any realistic choice of (3. A very steep slope (3 = 3.5 (shown by the dashed line) fits the 
observations, but would produce unrealistically large outflow velocities incompatible with 
the data. Moreover, we will show later that the cold shell becomes unstable when the slope 
of the density profile becomes steeper than (3 ~ 2: at that point fragments decouple from 
the hot wind. Thus, the steepest physical dependence of v w on the SFR is q ~ 1/3 when 
(3 = 2, while the best fit to the data requires q ~ 2/3. 

Guided by the aforementioned arguments, in this paper we propose that the different 
observed wind velocities correspond to different starburst evolutionary times, and are not 
directly related to variations of their mean SFR or other internal properties. The propor- 
tionality of the wind velocity and SFR is incidental. It reflects the fact that both quantities 
decrease as the time elapsed from the starburst increases. 

We have organized the paper as follows. In §|2]we describe our model analytically assum- 
ing a power law density profile for the gas. In § 01 we assume a more realistic density profile 
of the gas and we compare the numerical results to the observations. In § H] we discuss the 
implications of our model and provide an explanation for the lack of evidences for infalling 
gas in LBGs. 
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2 SCALING ARGUMENTS 

To understand some basic features of LBG outflows we first resort to scaling arguments. 
We closely follow the method outlined in Ferrara & Tolstoy (2000) and Madau, Ferrara & 
Rees (2001) to which we refer the interested reader for the details. Assume that ambient gas 
pressure, gravity, and cooling can be neglected; moreover, we take a constant mechanical 
luminosity for the starburst, L w (the case L w ^ const can be explored but there is no 
analytical solution). Then the outflowing (thin) shell radius evolution is described by: 

We further make the hypothesis that the halo gas density profile can be approximated (at 
least piece-wise) by a power-law profile: 
RY P __ ^0.22^ _ 



P9{r) = Pc {r~J =/ H~rJ x - p ° x > (2) 

with x = R/R V i r and p c is the central density. As for a Navarro, Frenk & White (1997, NFW) 
dark matter halo baryons have a core radius R c = (0.22/c)i?^ r , where c is the concentration 
parameter, we have used such relation in eq. 121 Note that po = Pc i n the case of a central 
core, i.e. if (3 = 0. By substituting p g into eq. [I] we find 



with q — 1/(5 — /?); we use the standard definition of R 
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We also take p c = 5 c {p g (z)) ~ 5 c VL b h 2 p crit {l + z) 3 . In the following, we adopt the value 
5 C = 2 x 10 4 , the appropriate overdensity of the core in a NFW profile. If we further define 1 
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eq. (j3J) can be casted in the very simple form: 

/ / \ 3q 

*.(*) = (-) =r 39 . (6) 



w 



The velocity of the shell is then given by 
at t 

1 We adopt the standard notation Y x = Y/10 x 

2 All lengths are normalized to R v i r , time to t w , velocities to R v i r /t w , masses to t^L/R 2 ^ 
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The dense shell is pushing on the rarefied interstellar medium (ISM) gas. The (external) shell 
interface becomes Rayleigh- Taylor (RT) unstable if the shell is accelerating. In this case a 
is parallel to v, and the vector goes from the heavy (shell) to the light (ISM) fluid and the 
RT instability develops. From MacLow & McCray (1988) simulations we note that the RT 
fingers are beyond the thin shell solution, i.e. the instability develops outside the bubble. 
The nonlinear evolution leads to the formation of clumps which detach from the shell; the 
hot gas wraps around them and escapes; at the same time, the ram pressure on the clumps 
decreases so they lag behind the shock. Hence, the onset of the RT instability is fixed by 
the condition that the acceleration of the shell becomes positive. The shell acceleration is 

^(r) = ^ = 3g(3g-l)r 3 - 2 . (8) 

The RT condition is satisfied when 3q(3q — 1) > 0; as q = 1/(5-/5) > for any reasonable 
density profile slope /3, this implies (3 > 2. Hence the instability starts to develop at the time 
Tj when the shell enters in the (3 > 2 part of the gas distribution. Before that it evolves in a 
(3 < 2 density profile (e.g. the baryonic core in a NFW profile or the inner part of a disk). 
The linear growth rate of the instability on spatial scale A = 27r//c is given by 



^rt(A,t) 



2- !)s (r) (A- l; 1 
A (A + l) 



(9) 



A 

the second equality descending by the fact that the density contrast between the shell and 
the ISM gas is A 3> 1. Note also that the instability grows first on the shortest scales and it 
diverges as A — > 0. At these small scales, the RT is stabilized by viscosity and/or magnetic 
fields. As we are interested in the shell fragmentation, we estimate the growth time on a 
scale A such that x c = X/Roi r , which by eq. ((HJ) is equal to rf q . The maximum growth rate 
occurs when the acceleration is largest and scale is smallest, i.e., at r = n. Then 

nV2 / , x/671- 1 (2-/3) | 1 

/6nq\(3q-l)\r-i= V , g K r~\ (10) 



r RT\Ji) = ^RTlTi 



2n\g s (Ti\ 



|5-/3| 



where < (3 < 2 is the value appropriate for the stable region. Given this range we find 
that the longest time corresponds to (3 = 0, for which T^T{ T i) — 0.8rj. So what is t(1 As 
already mentioned, the baryonic core in a NFW halo has size x c = 0.22/ c, hence for the 
(3 — (q — 1/5) solution and for c = 5 we find 

Ti = X V3? = x 5/3 = ( .22/c) 5/3 = 5.5 x 10~ 3 . (11) 

At the time -ry = Tj + trt = l-8xj, the instability becomes nonlinear and the shell fragments 
in dense clumps the move almost ballistically in the gravitational field of the dark matter 



9(x) = = _ M « (12) 
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halo from there on. Their precise mass/number is difficult to calculate, but to a first approxi- 
mation Mf ~ (it /6)p s f l 8R^ h , where 5R s h is the shell thickness calculated at 77; as now, we fix 
Mf = (J,M s h with \x < 1. The mass of the shell is approximately equal to the swept-up mass 
in the halo. To write the equation of motion for the fragments (note that this is independent 

of the mass of fragments) we first derive the nondimensional gravitational acceleration: 

2 

R 2 R vir \t ff J F(c)x 2 ' 

where t ff ~ [AnG p crit tt m (l + z) 3 (200/3)] _1/2 w 0.32[(1 + ^)/4]~ 3 / 2 Gyr is proportional to 
the halo free fall time. We have assumed a NFW profile with 

P{X) = cx(l + cx) 2 ' (13) 
5 = (200/3)c 3 /F(c), and 

F(f)= 1x1(1 + *)-^. (14) 
It follows that fragment velocity is governed by 

As cx <C 1, F(cx) « (c:r) 2 /(l + cx), then 3 

" ^ = _ 2 (*.y^/- _ A r (16) 

^(x/) \*/// -T (c) A/ (1 + CX) Jx f (l + CX) 

where Xf = x 8 (rf) = (1.8rj) 3<? , v s (x/) = 3gr^ 9_1 = 3g(1.8r,j) 3 ' ? ~ 1 and the constant A has 
value 

A = 2f— ) -^- = 2820. (17) 

Note that q ~ 1/5 (/? = 0) before r i; but it increases to g ~ 1/3 (/? = 2) after that, due 
to the steeper profile initiating the instability. We neglect that complication, by noting that 
the time-weighted mean is q = 0.3 ~ 1/3 when computing Xf and v s (rf). Thus, we obtain 
Xf « Tf and v s {rf) ~ 1. The resulting expression for the velocity is 

vj(x) = v s (x f ) 2 + A(x f - x). (18) 
Substituting the expression for v s (xf) above, and using g = 1/3 we find 



v f (x) = yjl + A(x f - x). (19) 



3 In the numerical results shown in Figurc|5]this approximation has been dropped. 
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The fragments are decelerated and they stop at xq = Xf + A^ 1 at a time tq obtained by 
solving 

x dx 



dr= I — ^ r = T f + 2A' 1 



l-Jl + A(x f -x) 



(20) 



evaluating this expression at x — xo, yields t = Tf + 2A~ 1 . For r > r , the fragments start to 
infall onto the galaxy unless, xq > 1, in which case they are lost to the gravitational potential 
and ejected in the IGM. Using the value of t w from eq. El , with [3 = 2 and c = 5 and the 
other fiducial values there, one obtains t w = 310L~ 40 Myr and A = 52(t w /tff) 2 = 2820L~ 2 / 3 . 
Then we find 

r = r f + 2A- 1 = 1.8r< + 7 x 10~ 4 L 2 J% = 0.01 + 7 x 10- 4 L^4 (21) 
x = Xf + A' 1 = Tf + A' 1 = 0.01 + 7 x 10" 4 L^4 ( 22) 

We conclude that mechanical luminosities L w ^ > 1.5 x 10 5 are required to eject the frag- 
ments (xq > 1) from a typical LBG galaxy. This conclusion is approximate as our formalism 
is valid in cx 1 (see eq. ITT)|) and it actually overestimates the velocity; nevertheless, it gives 
a first rough estimate of the escape of fragments. Are these mechanical luminosities consistent 
with the observed LBGs? The mechanical luminosity (stellar winds and SNe) for a continu- 
ous star formation rate of 1M yr" 1 , with M low = 1M , a = 2.35, M up = 100M Q , Z = 0.008, 
is logL = 41.8 (this slightly top-heavy IMF seems to be favored by the results obtained 
for high-redshift galaxies by Malhotra & Rhoads 2003). Assuming that a fraction rj rs 10% 
(Veilleuxei al. 2005) of this energy goes into kinetic energy, and the rest is radiated away, we 
find that L W ^ Q = 19(M*/30 M yr -1 ). Thus an enormous SFR of about 230 solar masses/yr 
would be needed (even more, given the overestimate of the velocity from eq. |16] ) to eject 
the fragment. The fact that cold, dense gas is very difficult to eject does not exclude that 
the hot gas can instead be expelled, as already noticed by Mac Low & Ferrara (1999). 

Suppose that after shell fragmentation the clumps move radially and keep constant size 
(this is a reasonable assumption as long as they are embedded in and pressure-confined by 
the hot bubble interior gas flowing into the halo). Then their covering factor of the source 

(which is independent of their mass) decreases with time as 

2 



C(r) 



Xf 


2 


"1.8V 


_x s (r)_ 




_x s (r)_ 



10-V- 2 , (23) 



that is, the covering factor has already dropped to 1% when x s = 10xf, which occurs at 
t = lOxf = 0.1, or in dimensional units, after 250L~ 1 / 3 Myr. 
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Finally, the hot shocked gas continues its expansion, leaving the clumps behind. Where 
does the shock stop? It stops when it reaches pressure balance with the IGM surrounding 
the halo, i.e. when its velocity has dropped to the IGM effective sound speed, which is in 
the formula below is assumed to be c s = 20 km s _1 . This happens approximately when the 
internal pressure L w t burst / R 3 is equal to (p g (z))c 2 s or, at z — 3, 

RT = [Lwtburst/{Pg(z))cl] ^ ~ ^ {j^j ^ (^) ^ M P C. (24) 

The previous argument neglects radiative losses, gravity and the possible infall ram-pressure; 
hence, such estimate represents an upper limit to R s (physical units). 



3 FULL MODEL 

The aim of this Section is to show quantitatively that the S03 data can be interpreted 
according to our model. By fitting the model to the data we attempt to infer some underlying 
properties of the LBG population. Hereafter, we adopt a more realistic description of the 
gas density profile with respect to the idealized case adopted in § 121 We also drop most of 
the simplifying approximations adopted previously. We solve numerically eq. (0) to calculate 
the evolution of the radius, velocity and acceleration of the shock front during the wind- 
driven phase. When the slope of the density profile becomes steeper than (3 = 2 the cold 
shell becomes unstable to RT instability and fragments. The subsequent evolution of the 
cold fragments is calculated by solving eq. (|T5|) and the time scale for shell fragmentation 
is calculated from eq. ©• The quantitative results found in this Section are robust as they 
are insensitive to the poorly known properties of the gaseous disk. We will show that our 
model is independent of the particular choice of the gas density profile, at least for "beta 
models" or exponential density profiles. The results are also insensitive to the values of 
physical parameters such as the disk scale height, core density etc., provided that they obey 
a given relationship. 

A gas in hydrostatic equilibrium in a NFW dark halo with virial radius r v ; r and concen- 
tration c, has a density profile that is well fitted by a beta model (Makino, Sasaki & Suto 
1998) 

n§ = [iTwkpF' (25) 

where, for a gas at the virial temperature, the core radius is R c ~ (0.22/c)R V i r and n c is 
the central gas density. It is likely that, due to gas cooling, the temperature of the gas is 
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lower than the virial temperature. In this case the core radius will be smaller and the central 
density higher. Independently of the assumed value of R c , if the gas mass is a cosmic fraction 
of the halo mass, M g = (Qi,/Qdm)M, by imposing mass conservation and disk geometry we 



We have assumed the following fiducial values for LBGs: total mass M = 2 x 10 11 M , 
concentration c = 5 (i.e.R c = 1.32 kpc). 

The data fitting is insensitive to the functional form of the density profile. We have found 
that adopting an exponential density profile n g = n c exp[— (R/R c ) a \ (where now R c is the 
scale height) the model fits the data equally well. The time scale for shell fragmentation 
depends on the shape of the profile. The evolution of the cold shell before fragmentation 
depends on the details of the profile but the subsequent evolution of the cold fragments 
is independent of the assumed profile. This happens because the fragmentation time-scale 
varies in such a way to keep the fragment velocity the same. We find that t^tI t i = 1.3 for 
the "beta model", trt/tj = 0.93 for the exponential profile with a = 1 and trt/tj = 0.66 
for a = 2. During the acceleration and fragmentation phases, the velocity of the shell 
increases only by a factor Av/vi of a few percent with respect to its initial value Vi at zero 
acceleration. This is consistent with eq. Q that can be recast in terms of the variables of 
interest: g s ~ Av/trt and x c ~ i^Tf. It follows that Av/vi ~ (27itrt / Vj) -1 <C 1. 

The velocity evolution of the cold fragments is degenerate for choices of the mechanical 
luminosity of the starburst L w that obey r/L w = Cn c R?J 2 , where C is a constant and 77 is 
the fraction of mechanical luminosity converted into kinetic energy of the cold shell (the 
rest being radiated away by the hot gas). Thus, further assuming eq. (|26[). it follows that 
the velocity of the fragments is proportional to r]L w R~ 1 ' 2 . The cold shell driven by the 
starburst is stable for only a few Myrs before it fragments and the clumps move ballistically. 
Thus the appropriate value of L w to use is that corresponding to this early phase, when the 
kinetic energy input produced by the starburst is dominated by winds from massive stars 
and is L w ^q = 1 — 5(M*/M yr -1 ), for stellar metallicities Z/Z Q = 0.4 — 1 (Leitherer et 
al. 1999). The value of rj, the fraction of mechanical luminosity converted into kinetic energy 
of the supershell, is not well known from observations. We adopt the fiducial value r\ = 10% 
(Vielleux et al. 2005). 

The LBG sample of S03 is grouped in four bins according to their Lya equivalent width. 



find 




(26) 
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Each bin contains about the same number of galaxies (about 200). In our interpretation, 
different bins describe a statistically steady-state time sequence and correspond to different 
times from the beginning of the starburst with typical time scale tburst- Thus, if we neglect 
differential selection bias effects in the data, we expect that equally time-spaced bins contain 
an equal number of galaxies of age: £j = {i/^)t burst f° r bins i = 1,2,3,4. The observed star 
formation rate of the galaxies in each of the four bins is well fitted by an exponential function 

M(t) = (29 M yr' 1 ) exp (-j^-) ■ (27) 

V ''burst ' 

The numerical value of tburst will be derived later by fitting the observations of the fragment 
velocities. It also follows from eq. (|27|) that the mean SFR integrated over the time scale of 
the burst is (M*) = 0.632 x M(0) M Q yr' 1 and the total stellar mass produced during the 
burst is M* = M(0)t burst « 10 9 M Q (t burst /30 Myr). 

In Figure we show the fit of our fiducial model (the "beta model" with R c = 1 kpc 
and n c — 0.1 cm~ 3 ) to the data from S03, shown by the symbols with errorbars. The panels 
from top to bottom show the time evolution of the star formation rate (fit by eq. [2*7j). the 
velocity of the hot gas (dashed line) and cold shell/fragments (solid line), the equivalent 
width of low ionization metal absorption lines, and the distance from the galaxy of the hot 
(dashed line) and cold fragments (solid line). The symbols with errorbars in the third panel 
refer to the normalized EW of Sill (circles), 01 (open circles), CII (squares) and Fell (open 
squares) calculated using eq. 

Good fits to the observed velocity and EWs of the metal lines are obtained for values of 
the starburst time scale 

tburst « (30 ± 5) Myr 

and for L w ^q = 1(M^/M yr _1 ) (appropriate for low stellar metallicity) we find 

where the second equality is derived by using eq. (J2T)|) that express mass conservation for 
the gas. As an example, if we instead assume an efficiency 77 ~ 20% we find R c ~ 100 pc, 
comparable to the disk scale height of the Milky Way. We also find that, independently of 
the assumed details of the gas density profile, the fragments start to fall back into the galaxy 
at tf a u « 60 Myr when they reach a distance from the disk of 2 kpc. Although on larger 
scales the outflow approaches spherical geometry, the fragmentation point is located close 
enough to the galactic center that the effects of the disk density stratification are important. 
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Figure 2. Outflow evolution as a function of time elapsed from the onset of LBG starburst. The symbols with errorbars show 
the S03 data (total of about 800 LBGs) grouped in four bins. The vertical lines show the time when shell fragments form and 
the hot gas freely escapes. From top to bottom the panels show the evolution of: (i) star formation rate, (ii) hot gas (dashed) 
and fragment (solid) velocities, (iii) normalized equivalent width of low ionization metal lines (Sill, OI, CII and Fell), (iv) 
radius of the hot bubble (dashed) and position of the fragments (solid) . 



4 FINAL CONSIDERATIONS 

We have already established that the cold clumps do not escape from LBGs. Therefore only 
two scenarios remain: either they disintegrate during the outflow or they eventually fall back 
onto the galaxy. In the second case the reason why we do not see infalling clumps is simple: 
we do not observe this evolutionary phase because it happens too late when the UV emission 
from the galaxy has faded and the galaxy cannot be found using the Lyman break selection 
criteria (selection techniques based on infrared colors might instead detect these galaxies 
long after the process has come to an end). More quantitatively, after about 50 — 100 Myr 
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the fragments will reach a maximum distance from the galaxy of about 2 kpc and than 
start falling back. According to our model the typical duration of the starburst phase is 
tburst ~ 30 ± 5 Myr. Although ages since the most recent episode of star formation are found 
to be < 40 Myr for 25% of the LBG sample (Shapley et al. 2001), the derived median value is 
320 Myr. A possible explananation for this tension could be that the stellar light is affected 
by the integrated star formation history of the galaxy (which could be more complex than 
assumed in population synthesis models used for the interpretation of the data), while the 
kinematic properties of the wind are primarily affected by the very last burst episode. To 
clarify this point additional study is required. 

During this time scale the SFR declines as exp [t/tb urst ] from 30 to 10 M yr _1 . During 
the starburst the fragments slow down from about 400 km s -1 to 100 km s _1 . The UV burst 
is not long enough to allow the observation of galaxies showing the infall phase of the gas 
fragments. In our model we have assumed that the fragments are stable and survive for at 
least ~ 30 Myr, that is approximately the age of LBGs grouped in the last bin. In other 
words the galaxies in the last bin show the low ionization metal absorption lines that arise 
from dense gas in the fragments. 

Let's now estimate the time scales for the destruction of the fragments. The two main 
processes that may disintegrate the fragments are shear instability (Kelvin-Helmholtz, here- 
after KH) and conductive evaporation. The time scale for the growth of the KH instability, 
neglecting saturation effects, is approximately 

tK « = 42 Myr fe) (20^) " iw) ^ ~ 50 - 100 Myr - < 28) 

where D = Pf/ph is the fragment-hot gas density contrast, Lf is the size of the fragments 
and Av their velocity relative to the flow of hot gas. We have estimated the fragment size 
assuming that is comparable to the thickness of the cold shell. When Av/c s > 0.6, the 
relevant case for our estimate, the instability growth rate saturates and eq. (|2*Hj) ceases to 
be valid (Vietri, Ferrara & Miniati 1997). In this case a good estimate of the minimum 
txH is derived by plugging the value Av = 0.2c s ~ 20 km s _1 in eq. ()28|) . At least for 
the largest fragments we find that txH — tfaii > hurst- Magnetic fields and cooling can 
increase the stability of the clumps, increasing our estimate of txH- A dimensional estimate 
of the evaporation time scale of the cold fragments embedded in the hot wind as a result of 
conductive heating gives 
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f N \ ( L \ f T \ ~^/ 2 
t mnd = 2^n / L? = 160Myr( 3xl0g f cm _ a ) M- (_) « 150-300 Myr, (29) 



.3 x 10 20 cm" 2 / \3 pc / UO 6 K, 
where Nf = UfLf rs (n g )xfR V i r /3 is the mean column density of the fragments (assumed 
to be equal to the column density of the cold shell). We have assumed Spitzer conductivity 
K e = 1.5 X 10 _7 T 5//2 , neglecting saturation effects and magnetic field (Spitzer 1953). 

The survival time to KH instability and thermal evaporation are about 100 Myr, a factor 
of three larger than the life-time of the star burst. The time scale for thermal evaporation 
is proportional to L 2 and for KH instability to Lf . We expect that the larger clumps may 
survive long enough to start falling back after tf a ii ~ 60 Myr but at this time the starburst 
has faded and the galaxy will not appear as a LBG with strong UV and Ly-alpha emission. 
It may be possible to test our model by selecting the faintest LBGs and searching for a 
subsample that has weak Lya emission. This subsample may show evidence of either low 
velocity outflows or gas infall. 

We would like to thank the referee, D. Weinberg, for insightful comments. 
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